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I.  INTRODUCTION 


Current  development  of  optical  amplifiers  has  advanced  to  the  stage  that  their 
widespread  use  in  lightwave  systems  is  certain  in  the  near  future  [1].  Receivers 
using  optical  amplifiers  have  been  shown  to  be  substantially  more  sensitive  than 
their  counterparts  [2].  Of  primary  importance  is  the  use  of  the  optical  amplifier  as 
a  preamplifier  in  a  direct  detection  receiver.  In  this  study,  a  unified  approach  to  the 
performance  evaluation  of  direct  detection  on-off  keying  (OOK)  lightwave  systems 
which  take  into  account  the  effect  of  postdetection  thermal  noise,  shot  noise,  and 
the  impact  of  nonzero  laser  linewidth  is  provided. 

The  direct  detection  OOK  receiver  to  be  analyzed  is  shown  in  Fig.  1.  The 
received  signal  is  amplified  by  an  optical  amplifier  of  gain  G.  The  optical  amplifier 
introduces  additive  white  Gaussian  noise  (AWGN)  with  zero  mean  and  power  spec¬ 
tral  density  (PSD)  N0/2  where  N0  =  Navhf(G  —  1)  [2,  12,  33].  The  parameter  Nsp 
represents  the  spontaneous  emission  factor  which  is  unity  for  an  ideal  amplifier,  h 
is  Planck’s  constant  (6.626  x  10-34  J.s),  and  /  is  the  frequency.  Because  the  optical 
amplifier  is  not  polarization  independent,  a  polarizer  is  needed  to  pass  the  desired 
signal  and  to  block  the  light  in  the  orthogonal  polarization.  In  this  study,  the  optical 
bandpass  filter  is  modeled  by  both  a  finite-time  bandpass  integrator  with  integration 
time  T'  [6,  9-11,  33]  and  a  Fabry- Perot  filter. 

For  the  finite-time  integrator,  the  noise  bandwidth  l/T'  of  the  bandpass  in¬ 
tegrator  is  chosen  to  be  the  same  as  the  noise  bandwidth  of  the  optical  filter.  The 
equivalent  lowpass  impulse  response  of  the  integrator  is  h{t)  =  (1  JT')[u(t)—u(t—T')\ 
where  u(t)  is  the  unit  step  function.  For  convenience  of  analysis,  it  is  assumed  that 
the  bit  duration  T  is  a  multiple  integer  of  T ',  that  is,  T  =  MT'  [33]. 
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Figure  1:  OOK  receiver  structure. 
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The  photodiode  has  a  responsivity  R  =  r}eq/hf  where  qe  <  1  is  the  quantum 
efficiency  and  q  is  the  electron  charge  (1.6  x  10-19  C).  The  output  current  of  the 
photodiode,  which  is  proportional  to  the  squared  envelope  of  the  receiver  signal,  is 
processed  by  a  detector  which  consists  of  a  low  noise  amplifier,  an  integrate-and- 
dump  filter  and  a  slicer  [33] . 

In  the  case  where  a  finite  time  integrator  is  employed,  the  amplifier  noise 
effect  is  most  pronounced  at  low  postdetection  thermal  noise.  With  a  higher  gain, 
the  amplifier  noise  is  dominant  and,  consequently,  the  net  amplifier  gain  is  reduced 
considerably.  In  the  case  where  a  Fabry- Perot  filter  is  employed  as  an  optical  filter 
at  the  receiver  model,  the  equivalent  lowpass  impulse  response  of  the  Fabry- Perot 
filter  can  be  well  approximated  by  an  RC  filter  for  both  single  ar  4  multi-channel 
within  the  frequency  range  \f  —  f0 1  <  FSR/20n  where  FSR  is  the  free  spectral 
range  and  f0  is  the  center  frequency  of  the  Fabry-Perot  filter.  For  example,  given 
FSR  =  3800  GHz,  the  approximation  works  very  well  for  | /  —  /0|  <  60.5  GHz.  In 
other  words,  the  effects  of  adjacent  channels  within  121  GHz  can  be  included.  The 
same  model  without  the  optical  amplifier  is  analyzed  in  [29]  and  it  is  seen  that  this 
model  enables  use  to  obtain  a  closed  form  analytical  expression  for  the  bit  error 
probability.  In  all  bit  error  probability  derivations,  whether  a  finite-time  integrator 
or  a  Fabry-Perot  filter  is  employed  as  the  optical  filter  in  the  receiver  model,  it  is 
assumed  that  all  channels  are  bit  synchronous  as  in  [29,  30]. 

This  thesis  is  organized  into  three  sections.  In  the  second  section,  the  receiver 
models  employing  finite-time  integrator,  and  the  analysis  is  a  reproduction  of  [33], 
Fabry-Perot  filter  for  single  channel  and  Fabry-Perot  filter  for  multi-channel  are 
analyzed.  Each  of  these  analyses  obtain  the  mathematical  framework,  bit  error 
probability,  and  the  numerical  results  deriving  the  detected  signal  envelope  and  its 
density  function.  The  last,  section  provides  the  conclusions. 
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II.  ANALYSIS 


The  analysis  is  organized  into  three  subsections.  All  three  provide  the  math¬ 
ematical  framework,  bit  error  probability  derived  from  the  statistics  of  the  white 
Gaussian  postdetection  thermal  noise  and  the  numerical  results  for  the  receiver 
model  considered.  The  receiver  model  to  be  analyzed  is  shown  in  Fig.  1. 

A.  FINITE-TIME  INTEGRATOR  MODEL 

The  optical  bandpass  filter  is  modeled  as  a  finite  time  bandpass  integrator 
with  integration  time  T' .  The  noise  bandwidth  1/T'  of  the  bandpass  integrator  is 
chosen  to  be  the  same  as  the  noise  bandwidth  of  the  optical  filter.  This  model  is 
used  because  it  is  analytically  tractable.  This  section  is  taken  from  [33]. 

1.  Mathematical  Framework 

For  mathematical  convenience  we  adopt  the  complex  envelope  notation 
of  a  real  OOK  signal.  Thus,  for  a  given  transmitted  bit  6,  the  corresponding  received 
signal  at  the  input  of  the  optical  bandpass  filter  in  Fig.  1  is  designated  as  follows: 

r,-(<)  =  y/GPbie}9^pr{t)  +  nc(t)  +  jn,(t)  (1) 

where  pr  =  l,0<t<T  and  zero  otherwise,  G  is  the  optical  amplifier  gain,  P  is  the 
peak  power  at  the  input  of  the  optical  amplifier,  9(t)  is  the  OOK  laser  phase  noise, 
and  nc(t)  and  n3(t)  are  the  independent  in-phase  and  quadrature  components  of 
the  additive  zero  mean  white  Gaussian  noise  representing  the  amplifier  spontaneous 
emission  noise.  The  PSD  of  nc(t)  and  nt(t)  is  N0/ 2  (we  use  No/2  instead  of  the  usual 
lowpass  PSD  N0  because  the  magnitude  of  the  signal  component  is  \JGP  bi  instead 
of  the  usual  y/2 GP  bi  for  complex  envelope).  Let  r  be  the  extinction  ratio  of  the 
laser  source  defined  as  the  ratio  of  the  transmitted  power  for  the  logical  zero  to  that 
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for  the  logical  one.  Then  6q  =  ^/r/(l  +  r)  for  the  logical  zero,  and  6j  =  +  r) 

for  the  logical  one. 

a.  Approximation  of  the  Detected  Signal  Envelope 
The  output  signal  of  the  optical  bandpass  filter  is  given  by 

si(<)  =  f;  /^Mt)  “  U(T  -  T'lr.fi  -  r)dr 

1  C7" 

=  pi  r'<‘-T><iT 

~  bSlrVi(T)iT  (2) 

The  photodetector  with  responsivity  R  detects  the  squared  envelope  f?|s,(<)|2.  This 
signal  current  plus  the  shot  noise  current  v,(f)  generated  by  the  photodiode,  and 
the  postdetection  thermal  noise  current  w(t)  are  scaled  by  1  /TRGP  and  integrated 
to  produce  the  decision  variable 

Yi  =  TRGP  L  R^t^dt  +  TRGpJ  Vi^dt 
+  TRGP  L  W^dt 


=  Xi  +  Vi  +  W 

(3) 

where 

*  =  trgpL  R^dt 

(4) 

Vi  =  TRGP  L  Vi(t)it 

(5) 

W  -  TRGP  L 

(6) 

We  model  w(t)  as  a  zero  mean  white  Gaussian  noise  current  of  spectral  density  Wo- 

It  is  seen  that  the  variance  of  of  u>(t)  is 

2  _  Wo 

°w  TR2G2P2 

(7) 
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The  shot  noise  current  v,(t)  can  be  modeled  as  a  zero  mean  Gaussian  noise  current 
with  spectral  density  Vo,  as  follows  [15]: 


The  variance  of  u,(f)  is  thus  given  by 


TRG2P 2 


£{M0l2} 


The  value  of  a\.  can  be  computed  from  (Al)  and  (A4)  as  follows: 

,  2 qJP$  (*PT  \  qMN^hfjG  -  1) 

v'  TRGP(jBT)2  \M  j  r2RG‘P‘ 


where  /3  is  the  laser  linewidth. 


The  statistics  of  Xt  are  very  difficult  to  obtain.  Therefore,  we  make 
the  following  approximation  of  Xt.  We  note  that  the  equivalent  lowpass  impulse 
response  h(t)  of  the  optical  bandpass  filter  is  assumed  to  be  time-limited  to  T’  = 
T(M.  Thus,  the  integrand  i2[si(fcT')!2,  k  =  1,2,  •  •  • ,  M  in  (4)  is  evaluated  at  t  —  kT' 
where  there  are  k  independent  and  identically  distributed  (iid)  samples  of  i?|s,(t)|2 
during  an  interval  T.  Thus  |si(fcT')|2  is  obtained  from  disjoint  integration  intervals 
where  the  phase  noise  9(t)  has  independent  increments  and  the  noise  nc(t)  and  n3(t) 
are  altogether  independent  over  disjoint  intervals.  Using  the  above  fact  we  can 
approximate  Xi  as  the  sum  of  M  iid  samples  T'R\si(kT')\2/TRGP.  That  is 


x-~mrP 

Substituting  (l)-(2)  into  (11)  we  obtain 


M 

E 

k=l 


where 


*«£  ~  +  Nck+jNsk 

k= i  vM 


Zfe  =  1  f1”  e^dt 
T  J(k-\)T> 
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(14) 


l  r 

Nck  =  - 7= —  /  nc(t)dt 

T'y/MGP  J(k-\)T<  cw 

1  tkV 

N,k  =  - I  n,(t)dt 

T'y/MGP  *v  ’ 


(15) 


The  random  variables  Nek  and  N,k  are  independent  Gaussian  random  variables  with 
zero  mean  and  variance  a1  given  by 


j2  N0  _Ntphf(G~  1) 

2TGP  2  TGP  1  ’ 

The  approximation  of  X{  in  (12)  is  equivalent  to  the  modehng  of  the  integrate-and- 
dump  filter  as  a  discrete-time  integration  that  sums  over  M  samples  taken  every  T' 
seconds  at  the  output  of  the  photodetector  [6,  13].  We  remark  that  both  X{  and  its 
approximation  have  the  same  mean  value.  If  js,(t)|2  has  a  constant  spectral  density, 
then  X{  and  its  approximation  also  have  the  same  variance.  This  happens  when  M 
is  large  and  the  spectral  density  of  |a,(t)|2  can  be  considered  constant  within  the 
bandwidth  1/T  of  the  integrate-and-dump  filter. 


b.  Probability  Density  Function  of  the  Detected  Signal 
Envelope 

From  (11)  we  observe  that  X{  is  approximated  by  the  sum  of  squares 
of  2 M  iid  Gaussian  random  variables.  Therefore,  X{  is  approximately  noncentral 
chi-square  distributed  with  the  following  conditional  probability  density  function 
(P<tf)  /*,(*.•  |7)  (18] 


fx/xih) 


TGP  ( MiA(M_1)/2  f  TGP(x?  +  tif-,/M)\ 
N.phf(G-l)[bh)  eXPl  N„hf(G- 1)  | 
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(17) 


(2  TGP^^~/M\ 

M~l[  Ntphf(G-l)  ) 

where  /ji/-i(-)  is  the  modified  Bessel  function  of  order  M  —  1,  and  7  is  the  value 
assumed  by  the  random  variable  defined  as  follows: 

r  =  Elz*|I  (is) 

*=1 

In  the  special  case  when  the  extinction  ratio  r  =  0,  then  60  =  0,  and  fx0(xo )  in  (17) 
is  a  chi-square  probability  density  function  (pdf)  given  by  [18] 

,  .  x  1  /  TGP  \M  w  .  f  TGPx 0  1 

fx‘M  ~  (M  -  iy.  \Nvkf(G -  1))  *°  exp{  N„hf(G-  1)1  (19) 

The  pdf  of  fxi(xih)  can  be  obtained  by  averaging  /x,(x;|7)  over  7.  Thus  knowing 
/r(7)  we  get 

fxiixi)  =  /  /^(ati|7)/r(7)d7  (20) 

Jo 

2.  Bit  Error  Probability 

For  a  detection  threshold  o,  the  bit  error  probability  conditional  on  the 
mean  x0  and  xx  assuming  a  combined  shot  noise  current  and  postdetection  thermal 
noise  current  spectral  density  Vo;  +  Wo  is  given  by 

ft(x0,x1)  =  ler/c(^)  +  ler/o(^)  (21) 

where  the  complimentary  error  function  erfc(-)  is  defined  as 

er/c(a)  =  ~  f  e~** dx  (22) 

y  T  Ja 

and  Oy.,  *  =0,1  are  the  variances  of  Y{  in  (3) 


aY. 


2  ,  2 
°V;  +  °W 

Ti?G2P2‘^{lS,(<)l  ^  TR2G 2 P2 


(23) 
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The  bit  error  probability  Pj  is  obtained  by  talcing  the  expectation  of  P(,(x0,  xj)  with 
respect  to  x0  and  Xj  by  using  (17) 

ft  = 

<24) 

The  optimal  threshold  that  minimizes  the  bit  error  probability  is  the  value  that 
satisfies  the  following  equation 

r erfc  (^)  f*°^dx° = r erfc  (^)  (xi  ^  <25) 

Prom  the  above  analysis  it  is  seen  that  the  evaluation  of  the  bit  error  probability  Pt, 
requires  the  knowledge  of  the  pdf  /r(7)  of  T  in  (18)  which  in  turn  requires  the  pdf 
of  \Zk\2  in  (13).  The  pdf  of  the  random  variable  |Z*|2  has  been  studied  extensively 
in  [13,  14],  [20,  21].  The  evaluation  of  P*  in  (24)  is  computer  intensive  even  in  the 
ideal  case  of  no  laser  phase  noise.  When  laser  phase  noise  is  taken  into  account,  the 
pdf  of  T  must  be  obtained  from  the  pdf  of  |Z* |2  via  an  M-fold  convolution. 

3.  Numerical  Results 

In  this  section  we  present  numerical  results  for  a  direct  detection  light¬ 
wave  system  with  the  following  parameters:  r  =  0.05,  1/T  =  500  Mb/s,  G  =  100 
and  1000,  and  WQ  =  10-24,  10-23,  10-22  (A2/Hz).  Figure  2  shows  the  bit  error 
probability  P*  versus  the  input  power  P  (dBW)  as  a  function  of  postdetection  noise 
spectral  density  Wo  for  a  receiver  without  the  filter  and  the  optical  amplifier.  Figure 
3  shows  Pj  versus  P  for  G  =  100,  M  =  10  without  phase  noise  and  with  phase  noise 
of  flT  =  1  as  a  function  of  WQ.  At  Pt,  =  10“15,  the  net  amplifier  gain  without  (with) 
phase  noise  is  19.6  dB  (19  dB)  for  W0  =  10‘22  A2/Hz,  18.6  (18  dB)  for  W0  =  10*23 
A2/ Hz,  and  16.3  dB  (15.7  dB)  for  Wo  =  10-24  A2/ Hz.  It  is  seen  that  the  amplifier 
noise  effect  is  most  pronounced  at  low  postdetection  noise. 
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Bit  Error  Probability 


Input  Power  (dBW) 


Figure  2:  Bit  error  probability  versus  input  power  for  a  system  without  optical 
amplifier  as  a  function  of  postdetection  noise  spectral  density  W0.  A:  Wa  =  10-22 
A2/Hz,  B:  W0  =  IQ"23  A2/Hz,  C:  W0  =  10'24  A2/Hz. 


Bit  Error  Probability 


Figure  3:  Bit  error  probability  versus  input  power  as  a  function  of  postdetection 
noise  spectral  density  W0  with  G  =  100,  M  =  10.  A:  W0  =  10"”  A2/Hz,  B: 
WQ  =  lO'23  A2/Hz,  C:  W0  =  10’24  A2/Hz. 


Figure  4  shows  Pj  versus  P  for  G  =  1000,  M  =  10  without  phase  noise 
as  a  function  of  Wo.  At  P*  =  10-15,  the  net  amplifier  gain  is  26.4  dB  at  W0  =  10-22 
A2/Hz,  22.1  dB  at  W0  =  10~23  A2/Hz,  and  17.2  dB  at  W0  =  lO'24  A2/Hz.  It  is 
seen  that  with  a  higher  gain,  the  amplifier  noise  is  dominant  and  consequently,  the 
net  amplifier  gain  is  reduced  considerably.  Comparing  the  two  systems  that  employ 
optical  amplifiers  with  gain  G  =  100,  and  G  =  1000,  respectively,  we  observe  that 
t  is  a  net  improvement  of  6.8  dB  at  W0  =  10~22  A2/ Hz,  3.5  dB  at  W0  =  10-23 
^  and  only  0.9  dB  at  Wo  =  10~24  A2/ Hz  for  G  =  1000  over  G  =  100.  Besides 
that,  Fig.  5  shows  the  optimized  threshold  versus  input  power  P  for  G  =  100, 
M  =  10,  30,  50,  and  100.  As  it  is  described  in  (25),  the  optimal  threshold  that 
minimizes  the  bit  error  probability  is  the  value  that  satisfies  the  equation  given 
in  (25).  As  we  can  see  from  Fig.  5,  optimized  threshold  value  decreases  as  input 
power  increases  for  all  M  values.  For  considered  values  of  M,  which  are  10,  30,  50, 
and  100,  the  optimized  threshold  appears  to  converge  to  about  0.5  for  higher  input 
power  values.  We  conclude  that  a  larger  gain  amplifier  should  be  used  when  the 
postdetection  noise  is  large,  and  a  smaller  gain  amplifier  should  be  used  when  the 
postdetection  noise  is  small. 

B.  FABRY-PEROT  FILTER  FOR  SINGLE  CHANNEL 

The  desired  signal  is  amplified  by  a  preamplifier  and  then  filtered  by  using  a 
Fabry-Perot  filter.  The  photodetector  has  a  responsivity  R  (A/W).  The  detected 
current  is  amplified  by  a  low  noise  amplifier  that  adds  a  postdetection  thermal 
noise  with  spectral  density  No  (A2/ Hz).  The  decision  variable  at  the  output  of  the 
integration  is  compared  to  the  threshold  in  order  to  be  able  to  determine  whether 
a  bit  zero  or  bit  one  is  present. 
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Bit  Error  Probability 


Figure  4:  Bit  error  probability  versus  input  power  as  a  function  of  postdetection 
noise  spectral  density  W0  with  G  —  1000,  M  =  10.  A:  W0  —  10~“  A2/Hz,  B: 
W0  =  lO"33  A2/Hz,  C:  W0  =  10~M  A2/Hz. 
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Normalized  Threshold 


Figure  5:  Optimized  threshold  versus  input  power  as  a  function  of  M  with  G  — 
1000,  postdetection  thermal  noise  W0  =  10~22  A2/Hz. 
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1.  Mathematical  Framework 

Again,  we  adopt  the  complex  envelope  (equivalent  lowpass)  notation  of 
a  real  00K  signal  for  mathematical  convenience.  For  single  channel  analysis,  we 
normally  ignore  the  intersymbol  interference  (1SI)  and  adjacent  channel  interference 
(ACI)  components.  Detailed  analysis  including  ISI  and  ACI  components  are  per¬ 
formed  in  the  multi-channel  case.  For  a  given  transmitted  bit  bt ,  the  corresponding 
received  signal  at  the  input  of  the  Fabry-Perot  filter  in  Fig.  1  is  designated  as 
follows: 

r,(t)  =  y/GPbofiPr{t)  +  nc{t)  +  jn»{t)  (26) 


where  G  is  the  optical  amplifier  gain,  P  is  the  peak  power  at  the  input  of  the 
optical  amplifier,  nc(t)  and  nt(t)  are  the  independent  in-phase  and  quadrature  com¬ 
ponents  of  the  zero  mean  AWGN  representing  the  amplifier  spontaneous  emission 
noise  (ASE)  and 

{1  0  <t<T 

(27) 

0  otherwise 

The  Fabry-Perot  filter  can  be  characterized  by  the  following  equivalent  lowpass 
transfer  function  [31,  32] 


»{f) 

H(f) 


1  —  p  1  —  A  —  p 

1  —  pe-i2*7/FS/i  1  —  p 


1  -  p 


1  —  p  cos 


-f  jp  sin 


1  —  A  —  p 

1  -  p 


(28) 


where  p  is  the  power  reflectivity,  A  is  the  power  absorption  loss  [A  =  0  for  an  ideal 
Fabry-Perot  filter)  and  FSR  is  the  free  spectral  range.  For  |/|  <  FSR/20ir  and 
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Thus,  if  the  signal  is  bandlimited  to  |/|  <  FSR/20n ,  we  can  actually  use  (29),  but 
with  the  frequency  covering  the  entire  frequency  spectral  range,  that  is,  we  can  truly 
approximate  (28)  with  a  single-pole  RC  filter  with  the  following  transfer  function 
and  impulse  response,  respectively, 


-  ,  \rf 

1+nr 

(31) 

h(t)  =  ce_cl  ,  t  >  0 

(32) 

The  magnitudes  of  H(f)  of  the  Fabry- Perot  filter  in  (28)  and  its  approximated  single¬ 
pole  RC  filter  in  (31)  for  p  =  0.99,  F  =  312.6,  B  —  12.16  GHz,  and  FSR  =  3800 
GHz,  remain  identical  and  attenuate  rapidly  while  the  phases  differ  as  the  frequency 
increases  [29]. 

The  above  approximation  is  valid  for  dense  wavelength  division  multi¬ 
plexing  (WDM)  analysis  when  the  filter  finesse  F  is  large  or  equivalently  the  FWHM 
bandwidth  B  is  small  since  the  equivalent  lowpass  signal  must  be  bandlimited  to 
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about  |/|  <  FSR/2Qjt.  When  |/|  >  FSR/20ir ,  the  magnitude  of  H(f)  is  very  small 
and  therefore,  the  effect  of  AC  I  beyond  this  frequency  range  is  negligible. 

a.  Derivation  of  the  Detected  Signal  Envelope 

Since  we  are  interested  in  the  detected  bit  6q,o  in  the  time  interval 


(0,T),  we  consider  the  output  preamplified  and  filtered  S,(t),  given  by 

Si(t)  =  SB{t)  0  <t<T  (33) 

where  Sg(i)  is  the  desired  signal,  which  is  the  convolution  of  the  desired  signal  and 
impulse  response  h(t)  in  (32),  as 

h(t)  =  ce~ct  t>  0  (34) 

5B(f)  =  VGP  bofi  [‘ h{t  -  r)dr  (35) 

SB(t)  =  v'GPb0,o( l-e-cl),  0  <t<T  (36) 


The  photodetector  with  responsivity  R  detects  the  squared  envelope  f?|5'l(<)|2.  This 
signal  plus  the  postdetection  thermal  noise  current  w(t)  are  scaled  by  1/TRGP  and 
integrated  to  produce  decision  variable, 

*  =  T&pC*S*tX*  +  TS5p£«,'>*  (37) 

Y,  =  Xi  +  W  (38) 


where 


*  =  trgpL  Wdt 

(39) 

w  =  trgpL  w(,)it 

(40) 

We  model  w(t)  as  a  zero  mean  AWGN  current  of  spectral  density  of  W0.  It  is  seen 


that  the  variance  of  o\j  of  w(t)  is 

2  W0T  W0 

~  T2(RGP)2  ~  T(RGP)2 


(41) 
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As  it  is  described  in  (37),  the  decision  variable  Yi  can  be  easily  found  for  zero  and 
one  bit  as  follows, 


Y0  =  X0  +  W  (42) 

Yx  =  Xi  +  W  (43) 

As  it  can  be  seen  from  above,  since  we  ignore  shot  noise  in  this  analysis,  the  variances 
of  Y\  and  Y?  are  equal  to  each  other. 

=  °yx  =  aw  (44) 


Xi  can  be  computed  from  (B.5)  and  found  as  (B.6) 

Xi  =  C  [i  -  gf(i  -  O  +  2^(1  -  «-")]  («) 

Using  the  fact  that  the  noise  nc{t)  and  n,(t)  are  altogether  independent  over  disjoint 
intervals  and  the  random  variables  Nc  and  N,  are  independent  Gaussian  random 
variables  with  zero  mean  and  variance  a\  computed  from  (B.18),  we  can  define 


variance  as 

2  N,phf(G  -  l)c 
u'  iMGP 


(46) 


b.  Derivation  of  Probability  Density  Function  of  the  Detected 
Signal  Envelope 

As  it  is  described  in  the  first  section,  Xi  is  approximately  noncentral 
chi-squared  distributed  with  the  following  probability  density  function  (pdf)  /Xj(x,) 
[18], 

f-M  -  R  (|)  (*f)  <«> 

where 


^  =  C  [l  -  J;(l  -  «'*■)  +  2^(1  -  «-JeT)]  (48) 
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*0.0  =  ^ 

l(  l-r)  fOT 

bit  one 

(49) 

Vo  =  ^ 

for 

(1  +  r) 

bit  zero 

(50) 

a?  «  MS-1  )c 
‘  4  MGP 

(51) 

In  the  special  case  where  the  extinction  ratio  r  =  0,  since  6q,o  =  0  for  bit  zero,  pdf 
of  fxi(xi)  in  (47)  becomes  a  chi-square  pdf  given  by  [18], 

/xo(lo)  =  (M  -  (52) 

2.  Bit  Error  Probability 

For  a  detection  threshold  a,  the  bit  error  probability  is  conditional  on 
the  mena  x0  and  Xi  and  given  by 

A(*„,  x.)  =  jer/c  +  ier/c  (^)  (53) 

where  er/c(  )  is  defined  as  in  (22)  and  cry.,  *  =  0, 1,  are  the  variances  of  Yi  in  (37) 

=  °Yt  =  40  =  aw  (54) 

We  can  represent  both  variances  as  cry,  since  they  are  equal  to  each  other.  The  bit 
error  probability  is  obtained  by  taking  the  expected  value  of  Pb(xo,  xi)  with  respect 
to  xo  and  X\  using  (47), 

Pb  =  4 1  er/c(vf^)/xo(Xo)<ilo+4 1  erfc{^~jfxAxi)dXi  (55) 

The  optimal  threshold  a  that  minimizes  the  bit  error  probability  is  the  value  that 
satisfies  the  following  equation  also  given  in  (25) 

f er/c  (^)  h'dx° = ferfc  Mz')iz'  (56) 


20 


3.  Numerical  Results 

In  this  section,  we  present  the  numerical  results  for  a  single  channel 
direct  detection  lightwave  system  employing  a  Fabry-Perot  filter  and  an  optical 
amplifier  at  the  receiver  model  in  Fig.  1  with  the  following  parameters:  bit  rate 
Rb  =  500  Mbps,  optical  amplifier  gain  G  =  100  and  1000,  photodiode  responsivity 
R  =  0.5,  FSR  =  3800  GHz,  c  =  38.4  GHz,  p  =  0.99,  thermal  noise  spectral  density 
Wo  =  10'22,  lO"23,  10-24  (A2/Hz),  and  M  =  10,  13,  20,  40  where  M  =  TJT ', 
T  is  the  bit  duration  (1/Rb)  and  T'  is  the  time  constant  of  the  optical  bandpass 
filter.  Since  bit  rate  Rb  —  500  Mbps,  bit  duration  T  is  2  ns.  If  we  chop  the  impulse 
response  of  the  filter  at  0.05  ns,  0.1  ns,  0.15  ns,  and  0.2  ns,  M  has  the  values  of 
40,  20,  13,  10,  respectively.  Figures  6,  7,  and  8  show  the  bit  error  probability  Pb 
versus  peak  input  power  P  (dBW)  as  a  function  of  M  for  the  thermal  noise  spectral 
density  W0  =  10“22,  10-23,  10"24  A2/Hz  with  G  =  100.  As  M  decreases,  it  is 
obvious  that  more  input  power  is  required  to  achieve  Pb  =  10~15  as  compared  to 
finite-time  integrator.  In  the  case  where  finite-time  integrator  is  employed  as  M 
decreases,  less  input  power  is  required  because  the  noise  variance  given  in  (16), 

■2  Ntphf(G-l) 

2  GPT 

is  independent  of  M  and  even  if  T'  is  varied,  it  has  constant  T.  Therefore,  varying  M 
does  not  effect  the  noise  variance  but  it  effects  the  signal  power.  For  the  Fabry-Perot 
filter,  the  variance  given  in  (46), 

2  N,phf(G-l)c 
4  MGP 

depends  on  the  M  value,  and  as  T'  is  varied,  M  varies  accordingly  (M  =  TIT'). 
So,  as  M  decreases,  noise  variance  gets  larger  and  more  power  is  required  to  achieve 
Pb  =  10~15.  As  it  can  be  easily  seen  from  these  figures,  the  role  of  T'  is  more 


Input  Power  (dBW) 


Figure  6:  Bit  error  probability  versus  input  power  as  a  function  of  M  with  G 
100,  WQ  =  10~22  A2/Hz  (Fabry-Perot  filter). 
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Figure  8:  Bit  error  probability  versus  input  power  as  a  function  of  M  with  G  = 
100,  W0  =  10’24  A2/Hz  (Fabry-Perot  filter). 
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significant  for  smaller  thermal  noise  values.  This  is  more  pronounced  with  a  high 
gain  amplifier  ( G  =  1000)  as  can  be  seen  in  Tables  1  and  2  and  Figs.  9,  10,  and  11. 
Detailed  comparison  is  as  follows: 

Table  1:  Comparison  between  M  =  10  and  M  =  40  with  G  =  100  (see  Figs.  6, 
7,  and  8). 


w„ 

(A’/Hz) 

M  =40 
(T  =  0.05  ns) 
(dBW) 

M=  10 
(T  =  0.2  ns) 
(dBW) 

Difference 

(dBW) 

10-22 

-71.0 

-70.0 

1.0 

10-23 

-75.2 

-73.0 

2.2 

io-24 

-78.9 

-73.8 

5.1 

Table  2:  Comparison  between  M  —  10  and  M  =  40  with  G  =  1000  (see  Figs.  9, 
10,  and  11). 


W0 

(A2/Hz) 

M  =  40 
(T  =  0.05  ns) 
(dBW) 

M  =  10 
{T  =  0.2  ns) 
(dBW) 

Difference 

(dBW) 

10-22 

-77.9  dB 

-73.8  dB 

4.1  dB 

»— * 
o 

1 

w 

-78.7  dB 

-73.9  dB 

4.8  dB 

Ci 

1 

O 

-78.8  dB 

-73.9  dB 

4.9  dB 
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Figure  10:  Bit  error  probability  versus  input  power  as  a  function  of  M  with 
G  =  1000,  W0  =  10~23  Al /Hz  (Fabry-Perot  filter). 
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When  we  compare  the  finite-time  integrator  and  Fabry-Perot  filter  with 
each  other,  Table  3  shows  that  the  finite-time  integrator  has  better  performance  than 
the  Fabry-Perot  filter  for  larger  T' .  This  improvement  of  the  finite-time  integrator 
over  the  Fabry-Perot  filter  gets  more  significant  with  high  gain  amplifiers  and  lower 
thermal  noise  values  as  displayed  in  Table  4.  Performance  margins  decrease  with 
increasing  thermal  noise  values. 

Figures  14  and  15  show  that  the  Fabry-Perot  filter  has  better  perfor¬ 
mance  than  finite-time  integrators  with  a  low  amplifier  gain  ( G  =  100)  for  low 
thermal  noise  values  (Wo  =  10-23  A2/Hz,  10-24  A2/Hz)  and  both  have  the  same 
performance  for  a  higher  thermal  noise  value  of  W0  =  10-22  A2/Hz  with  smaller  T', 
which  implies  a  higher  values  of  M  due  to  the  relation  of  M  and  T'  (Af  =  T/T')  for 
a  constant  bit  rate.  This  improvement  of  the  Fabry-Perot  filter  over  the  finite-time 
integrator  is  more  pronounced  with  a  high  gain  amplifier  ( G  =  1000)  as  shown  in 
Figs.  16  and  17  and  tabulated  in  Tables  5  and  6. 

As  a  result,  the  choice  of  T'  and  consequently,  the  value  of  M  plays 
a  great  role  for  the  systems  employing  the  Fabry-Perot  filter  and  the  high  gain 
amplifier  when  the  postdetection  thermal  noise  is  small.  The  systems  employing 
finite-time  integrators  and  high  gain  amplifiers  have  better  performance  than  those 
employing  the  Fabry-Perot  filter  when  T'  (time  constant  of  the  filter)  is  large  and 
thermal  noise  is  small.  The  systems  with  the  Fabry-Perot  filter  and  the  high  gain 
amplifier  have  better  performance  than  the  finite-time  integrator  when  T'  is  small 
and  thermal  noise  is  low. 
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Table  3:  Comparison  between  FTI  and  FP  for  G  =  100,  r  =  0,  and  M 
Figs.  6,  7,  8,  and  12). 


I 


Wo  (A2/Hz) 

FTI 

(Input  Power) 
(dBW) 

Fabry- Perot 
(Input  Power) 
(dBW) 

Remarks 
Improvement  of 
FTI 

io-22 

-71.1 

-70.0 

1.1 

io-23 

-75.3 

-73.0 

2.3 

10'24 

-78.5 

-73.8 

4.7 

Table  4:  Comparison  between  FTI  and  FP  for  G  =  1000,  r  =  0,  and  M  =  10  (see 
Figs.  9,  10,  11,  and  13). 


W0  (A2/Hz) 

FTI 

(Input  Power) 
(dBW) 

Fabry- Perot 
(Input  Power) 
(dBW) 

Remarks 
Improvement  of 
FTI 

io-22 

-78.5 

-73.8 

4.7 

io-23 

-79.6 

-73.9 

5.7 

10-24 

-79.8 

-73.9 

5.9 

Input  Power  (dBW) 


Figure  13:  Bit  error  probability  versus  input  power  as  a  function  of  postdetection 
thermal  noise  Wq  with  optical  amplifier  gain  G  =  1000,  extinction  ration  r  =  0, 
M  =  10.  A:  10-2  A2/Hz,  B:  10-23  A2/Hz,  C:  10-24  A2/Hz  (finite-time  integrator). 


Bit  Error  Probability 


Input  Power  (dBW) 


Figure  15:  Bit  error  probability  versus  input  power  as  a  function  of  postdetection 
thermal  noise  Wq  with  optical  amplifier  gain  G  =  100,  extinction  ration  r  =  0, 
M  =  50.  A:  10~"  A2/Hz,  B:  10’23  A2/Hz,  C:  10"24  A2/Hz  (Fabry-Perot  filter). 
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Figure  16:  Bit  error  probability  versus  input  power  as  a  function  of  postdetection 
thermal  noise  Wo  with  optical  amplifier  gain  G  =  1000,  extinction  ration  r  =  0, 
M  =  50.  A:  10"®  A2/Hz,  B:  10-23  A2/Hz,  C:  10-24  A2/Hz  (finite-time  integrator). 
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Bit  Error 


Input  Power  (dBW) 


Figure  17:  Bit  error  probability  versus  input  power  as  a  function  of  postdetection 
thermal  noise  Wq  with  optical  amplifier  gain  G  —  1000,  extinction  ration  r  =  0, 
M  =  50.  A:  10-22  A2/Hz,  B:  10’23  A2/Hz,  C:  10~24  A2/Hz  (Fabry-Perot  filter). 
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Table  5:  Comparison  between  FT  I  and  FP  for  G  —  100,  r  =  0,  and  M  —  50  (see 
Figs.  14  and  15). 


Wb  (A2/Hz) 

FTI 

(Input  Power) 
(dBW) 

Fabry- Perot 
(Input  Power) 
(dBW) 

Remarks 
Improvement  of 
Fabry-Perot  is 

w 

« 

1 

O 

-71.1 

-71.1 

0 

io-23 

-75.2 

-75.4 

0.2 

10-24 

-77.8 

-78.5 

0.7 

Table  6:  Comparison  between  FTI  and  FP  for  G  =  1000,  r  =  0,  and  M  —  50  (see 
Figs.  16  and  17). 


W0  (A2/Hz) 

FTI 

(Input  Power) 
(dBW) 

Fabry-Perot 
(Input  Power) 
(dBW) 

Remarks 
Improvement  of 
Fabry-Perot  is 

io-22 

-77.7 

-78.5 

0.8 

IQ-23 

-78.35 

-79.4 

1.05 

10-24 

-78.4 

-79.5 

1.1 
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C.  FABRY-PEROT  FILTER  FOR  MULTI-CHANNEL 

In  the  previous  section,  we  considered  only  the  desired  bit  signal  rather  than 
the  ISI  and  ACI  components  since  we  dealt  with  a  single  channel.  In  this  section, 
we  take  into  account  the  ISI  and  ACI  components  and  our  analysis  focuses  on  the 
multi-channel. 

1.  Mathematical  Framework 

For  convenience,  we  designate  channel  0  as  the  desired  channel  and  chan¬ 
nel  k  as  the  adjacent  channel  where  k  =  —M/2,  •  •  • ,  -1,1,-  •  • ,  M/2  and  M  is  an 
even  integer.  We  consider  the  equivalent  lowpass  (complex  envelope)  data  signal  in 
channel  0  and  channel  k  as  follows, 


o 


MO  =  KiPrit-iT) 

*=-6o 

(57) 

Mi)  =  £  ~  IT) 

(58) 

i——L 

where  T  is  the  bit  duration,  6o,.e{0, 1}  is  the  bit  in  channel  0  in  the  time  interval 
(iT,  (i-f  1)T),  6/t,/c{0,  eJ*}  is  the  £*h  bit  in  channel  k  in  the  time  interval  (£T,  (£+l)T), 
<t>ic  is  a  phase  offset  between  channel  0  and  channel  k,  which  is  assumed  to  be 
uniformly  distributed  in  (0,27t)  radians,  and  Wk  is  the  frequency  spacing  between 
channel  0  and  channel  k  in  radians. 

The  function  pr(t  —  iT)  is  defined  as 


r  i 


pr(t  -  iT)  =  < 


iT  <t  <  (i  +  l)T 


(59) 


l  0  otherwise 

The  non-negative  integers  L0  and  L  in  (57)  and  (58)  represent  the  number  of  bits 
in  channel  0  and  k,  respectively,  that  precede  the  detected  bits  6q,o-  The  received 
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WDM  equivalent  lowpass  signal  at  the  input  of  the  Fabry-Perot  filter  is  given  by, 

M/2 

r(t)  =  VGPbo(t)+  £  y/GPbk(t)  (60) 

kmU)  3 
*>*0 

where  P  is  the  received  optical  power  and  G  is  the  gain  of  the  preamplifier. 

As  we  discussed  in  the  previous  section,  the  impulse  response  of  the 
Fabry-Perot  filter  is  approximated  by  a  single-pole  RC  filter  of  which  the  impulse 
response  is  given  in  (32)  as, 

h{t)  =  cc~ct,  t>  0  (61) 

a.  Derivation  of  the  Detected  Signal  Envelope 

The  output  preamplified  and  filtered  signal  S,(<)  is  given  as, 

Si(t)  =  SB(t)  +  SISi(t)  +  SaciW  ,  0  <  t  <  T  (62) 

where  SB(t)  is  the  desired  bit  signal,  Sisi{t)  is  the  intersymbol  interference  signal, 
and  is  the  adjacent  channel  interference  signal.  These  signals  are  evaluated 

to  obtain  output  filtered  signals  by  using  (60)  and  (61), 

SB(t )  =  y/GPbofl  f*  h(t  —  T)dr 
Jo 

=  VGP  bovil  -  e-*) ,  0  <t<T  (63) 

Sisi(i)  =  VGP  ^2  &o,t  /  W  ~  T)dr 

i=Lo  J,T 

=  y/GPe~ct  2  bo,i  (e(*'+1)cr  -  eicT)  ,  0  <  t  <  T  (64) 

»=Lo 
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£ic/(0  =  Vgp 


m/2  rr  -i 

£  £*“/« 

km-M/2  ' 


h(t  -  t )e?WkT dr 


■¥biaj\(t-ry-'’dr^ 


M/2 

=  sfGPt £ 


km-U/2 


U<n,L  ' 


(c+jwk)(t+l)T  _  -(c+jw*)/T 


+6fc,o(e(c+^‘) 


-o}. 


0  <  t  <  T 


The  Fabry-Perot  filtered  output  Si(t)  is  detected  by  the  photode¬ 
tector.  The  detected  current  signal  is  amplified  by  a  low  noise  amplifier  (LAN)  that 
contributes  a  postdetection  thermal  noise  u>($)  with  spectral  density  N0  (A2/Hz). 
The  decision  variable  at  the  output  of  the  integration  is  compared  to  a  threshold  o 
to  determine  whether  bit  zero  or  bit  one  is  present.  The  decision  variable  Y  is  given 
by  (37)  as, 

Y  =  Xi  +  W  (66) 


where 


Xi  = 


W  = 


tSgp  £  R's‘Wdt 

hpJJw(t)dt 


TRGP 


We  obtain  the  signed  component  Xi  by  substituting  (62),  (63),  (64),  and  (65)  into 
(67)  as  a  function  of  cT  and  WkT,  which  represents  the  impact  of  IS  I  and  ACI, 
respectively.  Xi  can  be  obtained  from  (C.19)  for  the  worst  case  analysis  as  follows 

[29]  i 


Xi  = 


TRGP  [TRGP  (*L  [>  -  3*1  -  <-‘T>  +  SrU  -  e-)] 
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where 


+6-  y?(i  -  + **  [5=0  -  + •■">] 

+  [l*|2  +  g;(l  -  e-‘T)fle{6}(i_  -  Vo) 


60,0  =  0 
6_  =  1 
6=1 


M/2 

v  _ 

1  *w 

* 

km-M/2  1  + 

k*0 

(f)Jj 

« 

for  bit  zero 

6o,o  =  1 

for  bit  one, 

for  bit  zero 

0- 

I 

II 

0 

for  bit  one, 

for  bit  zero 

6  =  0 

for  bit  one 

(69a) 


in  the  worst  case  analysis. 

For  the  exact  case  analysis,  Xi,  which  is  derived  in  [29]  for  the  case 
without  the  optical  amplifier,  can  be  modified  for  our  case  as  follows, 


Xi 

Xi 


TRGpJo  R\SiW2dt 

TRGP  [TRGPblo  [*  “  &  (X  “  c  CT)  +  2cT 

+trGP-L  (l  - 


-1 


-  e 

t2 


E  Vi  (e,i+,)',■  -  e“T) 

Li.io 


+TRGP2tf  (*  - 


M/2  _i 


y'  y'  biej 


km-M/2 

k*0 


^e(cr+iu>*r)(<+i)c(cT+iu»*T)rj 
,  M/2  M/2 

+TRGP±  E  E 


bk,Obm,  0 


cT  ^  ^  (1  +jwkT/cT)(l-jwmT/cT) 

km-U/2  mm— M/2 


*  ¥0 


mft  0 
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{e-(cT-juikT)  _  i  e-(cT+jwmT)  _  i  1 

cT  +  1  -  jwkT/cT  +  1  +  jwnT/cT  +r1_ 


:-)]} 


i  (  M/2  M/2  _i 

+2(?/>T— Ke  j  E  E  E 


km-U/2  mm  —M/2 
0  mffO 


'L(l  +jwkT/cT)(\  —  jwmT/cT) 


(■ 


>(c7’+ju)*T)(#+l)  _  (cT+ju>tT)/ 


) 


1  —  g-(cT+jtumT)  i 

-K1- 


1  +  /cT  2 

+GPTVo4:  (l  +  e-icT  -  2£-eT)  E  Vi  (el'+,,‘r  -  eitT) 


r- )]} 


i=-Lo 


+GPTbofi^Rt^  (l  +  e~2cT  -  2e~cT) 

M/2  -i  . 

V  V  . °k'e  (  JcT\jwkT)(t+l)  _  JcT+jwkT)t\ 

£l1  +jw„T/cT'c  > 


km-M/2 

k*° 


+2  Y  V,o  l'e™T-l  T  1  1  T 

^  1  +  jwkT/cT\jwkT/cT^  ~  ~ 

km-M/2 
k*0 


2  2 


e-(cT-jwkT )  _  i 

1  —  jwkT/cT 

M/2  _! 


)} 


+  ^T  .E  V.,  (e^T  -  e'‘T) 

is-Lo 


M  £  lurfesfC--*’) 


km-M/2 

k*° 


^e(cT+jwkT)(t+ 1)  _  c(cT+iti>*T)^ 


M/2 


,  bkfi  ( 2(1  —  e 

r  ^  l+jwkT/cT\  1  -j 

km-M/2 


—  e-(eT~J,"*T)) 

'  +  e*2tT  -  1 


jwkT/cT 


)} 


JkflO 


(69b) 


Using  the  fact  that  the  noise  nc(t)  and  n,(t)  are  independent  over  disjoint  intervals 
and  the  random  variables  Ne  and  N,  are  independent  identically  distributed  (iid) 
Gaussian  random  variables  with  zero  mean  and  the  variance  af  computed  from 
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(B.18)  and  given  by  (46),  we  can  employ  the  variance  as, 


<r?  =  as  described  in  (46)  (70) 

4  MGP 

b.  Derivation  of  Probability  Density  Function 

As  it  is  explained  in  the  first  and  second  section  as  well,  A,  is 
noncentral  chi-square  distributed  with  the  following  pdf  also  given  in  (47) 

*■«>  =  (?)  M1/1  («) 
where  02  is  described  as  in  (69) 


II 

o 

£ 

^(1+r)  f°r 

bit  one 

(72) 

II 

o 

tljr— — r  for 

V  (Hr) 

bit  zero 

(73) 

*1  = 

#**/(<?“  1)« 
4  MGP 

as  in  (46) 

(74) 

The  special  case  where  the  extinction  ratio  r  =  0  introduces  the  chi-square  pdf  given 
by  (52)  and  [18], 


/* 'o(*o)  — 


TxP-'e-*0'20* 


~  (M  -  l)!(2og)"’ 0  v 

2.  Bit  Error  Probability 

For  a  detection  threshold  a  and  an  ISI/ACI  pattern  6  =  {6^.*,  hj}  where 
i  =  —Lq,  •  •  • ,  0  and  k  =  —M/2,  •  •  • ,  M/2(k  ^  0),  the  bit  error  probability  is  defined 
as  in  (53), 

n  m  _  1 _ (a  “  xo(*)N  ,  1 ,  f *l(&)  ~ 


ftW=rr/c(a^1)+Jer/c(i 


where 


2  2 


T(RGP)2 


Bit  error  probability  for  the  worst  case  analysis  and  exact  analysis  are  obtained  by 
taking  the  expected  value  of  Pt>(x 0,  Xj)  with  respect  to  x0  and  Xj  by  using  (69a)  and 
(69b),  respectively. 

Pt  for  the  worst  case  analysis: 

= \ f  ■* (^r)  '-<*»** + if  -*  (£f^)  (78) 

Pb  for  the  exact  analysis: 

h  -  h\r  /:  er/c 

b  \rrer/c  (xi(^i~q)  /-(•»»*« 

The  optimal  threshold  a  that  minimizes  the  bit  error  probability  is  the  value  that 
satisfies  the  following  equation, 

fer/c  (^^r)  /”(io)<'i° = r er/c  i^r)  Mx')dx'  (79) 

3.  Numerical  Results 

In  this  section,  we  present  the  numerical  results  for  the  system  employing 
the  optical  amplifier  and  the  Fabry- Perot  filter  in  the  receiver  model  for  multi¬ 
channel  with  the  following  parameters:  bit  rate  R*  =  2.56  Gbps,  optical  amplifier 
gain  G  =  1000,  free  spectral  range  FSR  =  3800  GHz,  M  =  3  and  4,  channel 
spacing  1  =  8,  12,  and  c  =  38.4  GHz.  In  our  analysis,  we  consider  that  the  signal  is 
band-limited  and  incorporate  the  degradation  caused  by  the  signals  in  the  nearest 
adjacent  four  channels.  Since  we  deal  with  multi-channel  in  this  section,  we  will  see 
the  effect  of  ISI  and  ACI  components  on  the  performance  of  the  system,  unlike  the 
single  channel  case. 
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For  the  worst  case  analysis  with  optimum  threshold,  bo, o  =  0  for  bit 
zero,  bo,0  —  1  for  bit  one.  Figure  18  shows  the  effect  of  ISI  and  ACI  for  M  =  3, 
cT  =  15,  and  I  =  8,  ’2  in  terms  of  bit  error  probability  and  input  power  (dBW).  As 
can  be  seen  from  the  graph,  the  single  channel  requires  —68.8  dBW  in  order  to  be 
able  to  achieve  bit  error  probability  Pj  =  10-15.  The  system  with  ISI  but  no  ACI, 
requires  —67.5  dBW  to  achieve  the  same  bit  error  probability,  which  corresponds 
to  a  1.3  dBW  power  penalty  as  compared  to  the  single  channel  case.  When  the 
adjacent  channels  are  taken  into  account  for  channel  spacing  I  =  12  and  8,  the 
power  penalties  axe  3.8  and  6.0  dBW,  respectively.  As  the  channel  spacing  between 
channels  gets  larger,  the  effect  of  ACI  becomes  less  significant. 

Figures  19  and  20  show  the  effect  of  the  value  of  M ,  which  is  the  ratio 
of  bit  duration  ( T )  to  the  time  constant  of  the  filter  impulse  response  ( T ').  As 
explained  in  section  2,  the  noise  variance  has  a  dependency  on  M  as  given  in  (46), 

_a  AW(G-l) 

1  4  MGP 

A  M  decreases,  the  noise  variance  gets  larger  and  more  noise  power  is  obtained 
accordingly.  As  can  be  seen  from  the  figures,  more  input  power  is  required  for 
M  =  3  rather  than  M  =  4  at  Pf,  =  10-15  as  expected  from  the  results  obtained  in 
the  single  channel  case  in  section  2. 

When  we  compare  the  worst  case  analysis  and  exact  analysis  with  each 
other,  as  we  expected,  exact  analysis  requires  slightly  less  power  than  worst  case 
analysis  in  order  to  obtain  the  same  probability  of  bit  error  P&  =  10" 15  as  shown  in 

Fig.  21. 
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Bit  Error  Probability 


G  =  1000 
M  »  3 


Figure  18:  Bit  error  probability  versus  input  power  for  Fabry-Perot  filter  with 
G  =  1000,  M  =  3,  W0  =  10-12  A2/Hz  (multichannel  worst  case  analysis). 
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Bit  Error  Probability 


Figure  19:  Bit  error  probability  versus  input  power  for  a  Fabry-Perot  filter  as  a 
function  of  M  with  G  =  1000,  W0  =  10“23  A2/Hz,  1  =  8. 
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Bit  Error  Probability 


Figure  20:  Bit  error  probability  versus  input  power  for  a  Fabry-Perot  filter  as  a 
function  of  M  with  G  =  1000,  W0  =  10“23  A2/Hz  and  I  =  12. 
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Bit  Error  Probability 


Figure  21:  Bit  error  probability  versus  input  power  for  Fabry-Perot  filter  with 
G  =  1000,  WQ  =  10-23  A2/Hz,  M  =  3,  I  =  8  (comparison  between  exact  and  worst 
case  analysis  for  multichannel). 
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III.  CONCLUSIONS 


We  have  presented  a  detailed  analysis  of  direct  detection  lightwave  systems 
employing  an  optical  preamplifier  at  the  receiver  and  derived  the  closed  form  ex¬ 
pression  for  the  bit  error  probability  of  WDM  systems  employing  on-off  (OOK)  as 
a  modulation  format.  We  have  considered  various  cases  in  which  the  receiver  model 
uses  either  a  finite-time  integrator  or  Fabry-Perot  filter  operating  in  a  single  chan¬ 
nel  or  multi-channel  environment.  We  have  taken  into  account  the  optical  amplifier 
noise,  the  postdetection  receiver  noise,  the  shot  noise,  and  the  effect  of  the  nonzero 
laser  linewidth. 

For  the  finite-time  integrator,  we  conclude  that  the  amplifier  noise  effect  is 
most  pronounced  at  low  postdetection  thermal  noise.  A  larger  gain  amplifier  should 
be  used  when  the  postdetection  thermal  noise  is  dominant,  while  a  smaller  gain 
amplifier  should  be  used  when  the  postdetection  thermal  noise  is  small  and  amplifier 
noise  dominates. 

The  performance  of  the  system  employing  the  Fabry-Perot  filter  in  a  single 
channel  environment  depends  on  the  value  of  M,  which  is  defined  as  the  ratio  of 
bit  duration  to  time  constant  of  the  impulse  response  of  the  filter.  When  we  chop 
the  filter  impulse  response  at  c  igh  value  of  Y,  which  requires  a  small  value  of  M 
based  on  the  relation  M  and  Y  for  a  constant  bit  rate  (M  =  T/Y),  we  conclude 
that  more  power  is  required  to  achieve  the  Pb  =  10-15  than  chopping  the  impulse 
response  at  a  lower  value  of  Y .  The  reason  for  this  is  that  the  noise  variance  of 
the  Fabry-Perot  filter  is  inversely  proportional  to  the  value  of  M.  As  M  decreases, 
more  power  penalty  is  paid  to  obtain  the  same  bit  error  probability. 
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When  we  compare  the  system  employing  a  finite-time  integrator  and  the  sys¬ 
tem  employing  a  Fabry-Perot  filter,  we  conclude  that  the  finite-time  integrator  has 
better  performance  with  a  high  gain  amplifier  under  low  thermal  noise  values  than 
the  Fabry-Perot  filter  when  we  chop  the  impulse  response  of  the  filter  at  a  high 
value  of  V .  For  a  larger  value  of  M,  we  conclude  that  the  Fabry-Perot  filter  has 
better  performance  with  a  high  gain  amplifier  and  low  thermal  noise  values  than 
the  finite-time  integrator. 

The  system  employing  an  optical  preamplifier  and  a  Fabry-Perot  filter  in  a 
multi-channel  environment  suffers  from  intersymbol  interference  and  adjacent  chan¬ 
nel  interference.  When  the  channel  spacing  is  large,  the  effect  of  ACI  is  small. 
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APPENDIX  B 


Derivation  of  X,  and  G, 


Xi  is  the  squared  envelope  signal  detected  by  photodetector  with  responsivity  R 


and  scaled  by  1/TRGP.  It  is  defined  as  in  (4) 

*  -  TMpJoR'S‘{t) t*  (B1) 

S,(t)  =  Sait)  =  y/GPbofi  f  h{t  -  r)dr  (B.2) 

h(t)  =  ce~ci ,  t>0  (B.3) 

Si(t)  =  VGPbofl  (l  -  e~ct)  ,  0  <  t  <  T  (B.4) 

*  =  tMp£R  (B.5) 

X  =  Vo  (l  -  |  (l  -  e-‘T)  +  ^  (1  -  e'**))  (B.6) 


It  is  described  in  the  first  section,  the  noise  nc(t )  and  n,(t)  are  independent  over 
disjoint  intervals  and  the  random  variables  Nek  and  Ntk  are  independent  identically 
distributed  Gaussian  random  variables  with  zero  mean.  Their  variances  in  (46)  can 
be  computed  as  follows, 


KW  =  7WpLh{t'T)n' (T)iT 

h(t)  =  hoo(<)PT'(0 

ii  o  <t<r 

0  otherwise 


Nc(t) 


7W?Lh~{t-T)nc{T)Mt'T)dT 

7M?Lh~{t-T)n‘{T)dT 


(B.7) 

(B.8) 

(B.9) 
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«c(r)  =  — 


i  rr  h 

y/MGPJ(n-i)r  00 

No 

2 

No  rr 


using  t'  =  nT'  —  r, 


2  2 
aNc  =  aN, 


2  _  _2 
^Nc  - 


2MGP  Jt< 

__No_[T> 
2MGP  Jo 

No 

2  MGP 


£kiw-ds) 

£  hUr)dr  <^pf  hlW 


(B.10) 


(B.ll) 


"  <  =  Mb£-«rh^T'-^  <B12) 


(B.13) 


where 


£-  fT'hl(T)di 

Jo 


(B.14) 


If  we  use  the  approximated  RC  filter  of  which  impulse  response  h(t)  given  as 


in  (B.3) 


where 


2  2 
°nc  =  °n.  = 


2  MGP 


(B.15) 


e  =  ['him 

Jo 

=  f\<x-«Yi 


£  =  J<l-e-!'r) 


(B.16) 


We  can  approximate  £  as  c/2  since  the  exponential  term  goes  to  zero  due  to  the 
high  value  of  2 cT'. 

£*>  *  f  (B.17) 
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Finally,  the  variable  is  found  as 


aNe  =  °N,  « 


2  2 

=  Of*. 


No  c 
2MGP  '  2 
N,phf(G  -  l)c 
4MGP 


(B.18) 


We  also  know  that  this  variance  value  can  be  represented  by 


aN,  =  in  Pd!  I54) 


(B.19) 


where  Ngp  (the  spontaneous  emission  factor)  is  unity  and 


FSR(  1  -  p) 


(B.20) 
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APPENDIX  C 


Data  signal  in  channel  k: 

MO  =  E  -  IT) 

(C.l) 

t=-L 

Data  signal  in  channel  0: 

M0-  E  k,PT(t-iT) 

(C.2) 

i=—Lo 

Bit  in  channel  k  in  time  interval  (IT,  (t  -f  1  )T) 

h,e  €  {0, 

Bit  in  channel  <f>  in  zth  time  interval  (iT,  (i  -f  1  )T) 

&o,»  €  {0,1} 


Detected  bit  in  (0,  T) 

bo,o  €  {0, 1} 

depending  on  where  bit  zero  or  bit  one  is  present. 

Desired  channel  =  channel  0 
Adjacent  channel  =  channel  k,  where 

M  M 

k  =  —  ,—1,1,  — even  integer 

Wk  —  FVequency  spacing  between  channel  k  and  channel  0 

2  nkl 
=  ~7F~ 


(C.3) 


< pic  =  Phase  offset  between  channel  k  and  channel  0  uniformly  distributed  over  (0,  2tt) 

1  0  <  t  <  T 


Pr(t)  = 


pr(t-iT)  = 


0  otherwise 

1  iT<t<  (i  +  1)T 
0  otherwise 


(C.4) 
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The  received  signal  at  the  output  of  the  preamplifier  and  at  the  input  of  the  Fabry- 
Perot  filter  (for  channel  <f>) 

M/2 

r(t)  =  VGPb0(t )+  £  VGPbi(t)  (C.5) 

km-Ufi 

k*0 

where  P  is  the  optical  power  and  G  is  the  optical  amplifier  gain. 

The  output  of  the  Fabry-Perot  filter  is 

r0(*)o  /  h{t  -  T)r(r)dr  (C.6) 

OO 

where  h(t)  is  the  equivalent  lowpass  impulse  response  of  the  Fabry-Perot  filter  of 
channel  0, 

h(t)  =  ce~ct  t  >  0  (C.7) 

,oo  m/2  f00 

r0(<)  =  VGP  I  h(t  -  r)fe0(r)dr  +  VGP  /  M*  “  T)bk(r)dT 

J  —  OO  J— OO 

km-M/2 
k*  0 


=  VGPbo,o  f°°  h(t  -  r)pr(r)dr 

J  — OO 

+VGP  £  &o ,i  /  h(t  -  t)pt(t  -  iT)dr 

■  r  J  —  OO 


i=~I»o 


M/2  0 


+'/GP  £  t,  /“  -  ’  )e>“*'pr(T  -  ffjir 


(C.8) 


km-M/2 


Since  we  are  interested  only  in  the  detection  interval  0  <  t  <  T,  we  need  to  evaluate 


S(t)  =  r0(t),  0  <t<T 


_ _  ft  - -  ~i  r{i+l)T  - -  '"1* 

S(t)  =  y/GPbo,o  h(t  -  T)dT  +  y/GP  £  K  .  h{t-T)dr  +  VGP  £ 

Jo  ,=-L0  JiT  ,  .. 


km-M/2 


(  £  r"T  Ht~  r)emTdr  +  b„,0  f  h(t  -  r )e*»MT 

\  [/=-L  J  JO 


(C.9) 
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S{t)  =  Ss{t)  +  Sici{t)  +  SACi{t) 


where 


SB(t)  =  VGP  boi0  f  h(t  -  T)dT 
Jo 

Sisi{t)  =  V GP  Y  bo.i  [  h(t  -  T)di 

«=-Lo  JiT 


M/2 


SACi(t)  =  VGP  Y 


kmM/2 

*yo 


r(‘+i)T 

E  /„  M<  -  r)e>*»Vr 

<=-L  JCT 


+b,:flJ\{t-Ty^rdT^ 

The  evaluation  of  the  integrals  are  as  follows 


Ss(t) 
Sisi(t ) 


M/2 


SAci(t )  =  y/GPcc~ct  Y 


1 


km-M/2 
k*  0 


C  +  jwk 


Y  bk,t  (e(c+iu'*)(<+1)r  -  e{c+jWk)T)  I  +  bkt0  (e{c+jWk)t  -  1 
t=-L  '  J  \  '  \ 


(C.10) 


(C.ll) 

(C.12) 


(C.13) 


VGPboflil-e-*),  0  <t  <T  (C.14) 

VGP  Y  V.  (e(i+1)cT  ~  eicT)  ,  0  <  t  <  T  (C.15) 


(C.16) 


Worst  Cast  Analysis 
bo,i  =  L0  =  oo 

S?si(t)  =  VGPb.e'*  0<t<T 


(C.17) 
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bk,t  =  bk,o  =  b,  L  =  oo 


_ _  W/J  eiw>‘t 

S?c,(t)  =  VGPb  £  - JS7  0<t<T  (C.18) 

km-U/J  1  + 

*»<0 

with  /  =  integer  >  0  and  in*  =  2irkI/T, 

X  =R  F  |5,(t)|2dt  (C.19) 

Jo 

By  using  (C.14),  (C.17),  and  (C.18),  the  evaluation  of  integrals  gives  the  value  of 
X  as 

*  =  5^?  +  5^0-.- <*)] 

+£  [j^(l  -  e'^)]  +  M-  [g=(l  -  +  e-™)\ 

+  [|6P  +  |=(1  -  e-'T)Jte{b}(b.  -  V,o) 


**0 


(C.20) 
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